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ABSTRACT 

To study the atomic, molecular and ionized emission of Giant Molecular Clouds (GMCs) in the Milky Way, we have initiated a Large 
Program with the Karl G. Jansky Very Large Array (VLA): ‘THOR - The H i, OH, Recombination Line survey of the Milky Way’. We 
map the 21cm Hi line, 4 OH lines, up to 19 Ha recombination lines and the continuum from 1 to 2 GHz of a significant fraction of 
the Milky Way (I = 15° - 67°, |J>| < 1°) at an angular resolution of ~ 20". Starting in 2012, we mapped 4 square degrees of the GMC 
associated with the W43 star-formation complex as a pilot study. The rest of the THOR survey area was observed during 2013 and 
2014. In this paper, we focus on the Hi emission from the W43 GMC complex. Classically, the Hi 21cm line is treated as optically thin 
with properties such as the column density calculated under this assumption. This approach might give reasonable results for regions 
of low-mass star-formation, however, it is not sufficient to describe giant molecular clouds. We analyzed strong continuum sources to 
measure the optical depth along the line of sight, and thus correct the Hi 21cm emission for optical depth effects and weak diffuse 
continuum emission. Hence, we are able to measure the H i mass of this region more accurately and our analysis reveals a lower limit 
for the Hi mass of M = 6.6_i g x 10^ Mq (vlsr = 60 - 120kms“'), which is a factor of 2.4 larger than the mass estimated with the 
assumption of optically thin emission. The Hi column densities are as high as Ah, ~ 150Mopc“^ x: 1.9 x lO^^cm^^, which is an 
order of magnitude higher than for low mass star formation regions. This result challenges theoretical models that predict a threshold 
for the Hi column density of ~ lOMepc^^, at which the formation of molecular hydrogen should set in. By assuming an elliptical 
layered structure for W43, we estimate the particle density profile. For the atomic gas particle density, we find a linear decrease toward 
the center of W43 with values decreasing from whi = 20cm“^ near the cloud edge to almost 0cm“^ at its center. On the other hand, 
the molecular hydrogen, traced via dust observations with the Herschel Space Observatory, shows an exponential increase toward the 
center with densities increasing to iih, > 200cm“^, averaged over a region of ~10pc. While at the cloud edge atomic and molecular 
hydrogen are well mixed, the center of the cloud is strongly dominated by H 2 emission. We do not identify a sharp transition between 
hydrogen in atomic and molecular form. Our results are an important characterization of the atomic to molecular hydrogen transition 
in an extreme environment and challenges current theoretical models. 
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1. Introduction 


It is well known that stars form in Giant Molecular Clouds 
(GMCs, [Mac Low & Klessen][2004| [McKee & Ostriker|[2007l 
Dobbs et al.|2014[|Offner et al. 2014|l, which primarily consist of 
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molecular hydrogen. Yet it is still under debate whether molecu¬ 
lar hydrogen is actually necessary for star formation or whether 
molecular hydrogen and stars form only under the same condi¬ 
tions side by side ( |Glover & Clark|[2012| l. The density within 
these clouds is high enough (particle density; n > 1 - 100 cm“^, 
column density: N > 1 - lOOMopc^^ ~ 10^® - lO^^cm^^) 
for molecular hydrogen to become self-shielded from the 
interstellar radiation field, which would otherwise dissociate 
the H 2 molecules (e.g. Dobbs et al.| 2014| l. Hence, molecular 
clouds form in the interior of large clouds of neutral hydrogen, 
which themselves are the environment of molecular clouds. 
Another open question is the fraction of the neutral hydrogen 
within molecular clouds and how this affects the physical 


conditions ( |Rrco & Goldsmith] |2010[ [Goldsmith & Li|[2005| l. 
Furthermore, the correlation between the transition between 
these two fundamental states of hydrogen and the corresponding 
physical conditions are also still unclear. 

As cold molecular hydrogen is challenging to observe directly, 
it is difficult to study its distribution in detail. Classically, the 
low-J rotational transitions of CO are used as a tracer for H 2 . 
However, recent simulations and observations show, that a large 
amount (~ 40%) of H 2 is CO-dark and therefore not well traced 
by CO ( Pineda et al.|2013] Smith et al.|2014| l. Another approach 
to study molecular hydrogen is via observations of thermal 
dust emission or dust extinction ([Lada et al. 2007t Molinarij 


jet al.|2010[|Kainulainen & Tan|20i3HKainulainen et al.|2013] l. 

On the other hand, the 21cm spin-flip transition of hydrogen 
offers a well-known method to measure the atomic gas content. 
Even though the 21cm line is well studied (e.g. Radhakrishnan 


et al.|197T} [Gibson et al.|2000[ [Taylor et al.|2003HGibson et al 


2005a[ Heiles & Troland 2003a|b| jStrasser & TaylOT 2004 
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disentangle the different contributions of the cold and warm 
Hi in emission and absorption for different spin temperatures 
and optical depths. On top of that, radio continuum emission 
at the frequency of the Hi emission line can also suppress the 
intensity of the observed H i emission, an effect that is especially 
significant for the Galactic Plane. 

To address cloud formation. Hi to H 2 formation as well as 
many other issues in ISM studies of the Milky Way, we initiated 
the THOR project: The Hi, OH, Recombination Line Survey 
of the Milky Way (Beuther et al. in prep.). We are using the 
Karl G. Jansky Very Large Array (VLA) in a large program to 
observe the Hi line at 21cm, four OH lines at 1612, 1665, 1667 
and 1720MHz, 19 Ha radio recombination lines (RRL) and 
the continuum from 1 - 2 GHz, which is split in 8 sub bands. 
These observations result in an angular resolution of 15 - 20". 
We have been awarded more than 200 hours observing time to 
map a significant fraction of the Milky Way (Galactic longitude 
/ = 15 - 67°, Galactic latitude \b\ < 1°) for the semester 2013A 
and 2014B. To test our observing strategy and data reduction, 
we began with a pilot study of a 2° x 2° held around W43, which 
was completed in 10 hours. The following paper describes this 
pilot study, which was observed in the semester 2012A. As the 
wealth of this data set is immense, we will focus on the H i line 
in this paper, but consecutive papers will focus on other aspects 
of the pilot study (e.g. Walsh et al. in prep., Johnston et al. in 
prep.). The full survey will be presented in Beuther et al. (in 
prep.). 

The held chosen for the pilot study is around the massive star 
forming complex W43 {I - 29.2 - 31.5°, \b\ < 1°). This region 
is situated at the intersection of the Galactic bar and the hrst 
spiral arm (Scutum-Centaurus Galactic arm, [Nguyen Luong| 
et al.|20lT} Carlhoff et al.|2013]l, leading to complex kinematic 


structures and possibly high star formation activity. W43 is 
referred to as a Galactic mini-starburst region ( [Motte et al.| 
2003 1 Bally et 'aL||2010 [) and shows a s tar-formation rate of 
0.1 Mq yr~‘ ( Nguyen I^uong et al.|201 1 1 or 5-10% of the star 


formation rate in the entire Milky Way. [Motte et al. ( [2014[ l found 
velocity gradients in CO and Hi 21cm emission. These velocity 
gradients could indicate large scale velocity Hows, which could 
cause the vast star-formation activity. [Motte et al.[ ( [2014| l also 
discussed that they do not hnd a threshold for the Hi column 
density, which is proposed by theoretical models ( [Krumhol^ 
|et al.|200^[2009| l. They argue, that we see several transition lay¬ 
ers of HI and H 2 along the line of sight and that such models are 
not suited to describe complicated molecular cloud complexes 
such as W43. The center of W43 harbors a large Hii region 


which is fueled by a Wolf-Rayet and OB star cluster (Liszt et al. 
[1993[ [Lester et al.|1985| [Blum et al.|1999| l. Furthermore, W43 
exhibits several high-mass starless molecular clumps, which 
are still in an early stage of star-formation ( Beuther et al.|2012t 
[Louvet et al.|2014) . Some massive dense clumps can potentially 
form Young Massive Clusters, progenitors of globular clusters 
( Louvet et al.|2014] Nguyen-Lu’o’ng et al.|2013] l. The BeSSeL 
survey (Bar and Spiral Structure Legacy Survey, [Brunthaler| 
[et al. [2017] [Reid et al.|2014] l has determined the distance to W43 
to be 5.5+0.4kpc from parallax measurements of methanol and 
water masers ( [Zhang et al.|[2014| ). This result has to be treated 
cautiously, as none of the four masers used for the parallax 
measurements are spatially directly associated with W43-Main. 



U [kA] 


Fig. 1 . UV coverage of one pointing (centered at RA 18:47:54.125746, 
Dec -03:28:42.90941, J2000) after 4x2 min observing time. 


2. Observations and data reduction 

2.1. VLA Observations 

We mapped a 2°X2° held around W43 (/ = 29.2-31.5°, |h| < 1°) 
during the 2012A semester with the Karl G. Jansky Very Large 
Array (VLA) in New Mexico in C-conhguration (Project 12A- 
161). As we used Nyquist sampling at 1.42 GHz with a primary 
beam size of 32', we needed 59 pointings to cover the 4 square 
degree mosaic. We chose a hexagonal geometry for the mosaic, 
similar to the VGPS survey ( [Stil et aL||2006| l. This results in a 
smooth areal sensitivity function with huctuations of less than 
1% in the inner region and a decreasing sensitivity toward the 
edges of our held. Each pointing was observed 4x2 min, which 
results in an overall observation time of 10 hours, including ~2 
hours overhead for hux, bandpass and complex gain calibration. 
The resulting uv-coverage for one pointing after 4x2 min ob¬ 
serving time is shown in Fig. The observations were done in 
two blocks each with 5 hours observing time in April 2012. We 
chose the quasar 3C286 as a Hux and bandpass calibrator and 
the quasar J1822-0938 as a complex gain calibrator, which was 
observed every ~ 13 min. Observing at L-band and using the new 
WIDAR correlator, we were able to simultaneously observe the 
Hi 21cm line, 4 OH lines (1612, 1665, 1667 and 1720 MHz) 
and 12 Ha RRL. For the pilot held, we had the spectral capa¬ 
bility to observe 12 Ha RRL, however, for the full THOR sur¬ 
vey we are able to observe 19 Ha RRL. The continuum, con¬ 
sisting of eight spectral windows between 1 and 2 GHz, was 
observed in full polarization. For the Hi 21cm line we used a 
bandwidth of 2 MHz with a channel width of 1.953 kHz. This 
results in a velocity range of +200 km s * and a channel spacing 
of Av 0.41 kms^k 

2.2. Calibration 

To edit and calibrate the data, we use CASA (version 4.1.0) with 
a modihed VLA pipelin^ (version 1.2.0). The pipeline does au¬ 
tomatic bagging for e.g. zeros or shadowing of antennas. Addi¬ 
tional bagging for Radio Frequency Interference (RFI) and bad 
antennas was done manually by hand. The pipeline also applies 
the bandpass, hux and complex gain calibrator to the data. We 
do not use Hanning smoothing and do not recalculate the data 
weights (’statwt’) at the end of the pipeline. The implemented 

* https://science.nrao.edu/facilities/vla/data-processing/pipeline 


Article number, page 2 of 


17 















































































S. Bihr et al.: THOR - The H i, OH, Recombination Line Survey of the Milky Way 


modifications help to improve the quality check, and we do some 
flagging and editing by hand subsequently. A full description of 
our quality check method will be presented in our forthcoming 
overview paper (Beuther et al., in prep). 


2.3. Imaging and deconvolution 


The HI emission and absorption covers a large range of angu¬ 
lar scales, which is challenging for the data reduction, as the 
VLA C-Array filters out most of the large scale structure. There¬ 
fore we combined our THOR data with the VGPS data (VLA 
D-Array and Green Bank Telescope data, |Stil et alT 2006| l, to 
overcome missing flux problems and to reconstruct the large 
scale structure. We tested three different methods to combine 
the THOR and VGPS data. First, we tried to combine the vis¬ 
ibility of the VGPS and THOR data and clean them together. 
Second, we tried the ’feather’ command in CASA and third, we 
used the images of the VGPS survey as a starting model (param¬ 
eter ’modelimage’) for the clean process in CASA. After testing 
these three methods, we choose the third method, as this method 
provides the best results considering noise, side lobes and re¬ 
covery of large scale structure. We compared point source peak 
fluxes of the combined images to the VGPS data to check for 
consistency. The difference is on the 10% level. The clean pro¬ 
cess was stopped at the 5cr noise level. 

We smoothed over three channels to reduce the noise, resulting 
in a velocity resolution of 1.24 km s“'. This was the best com¬ 
promise between computational time, noise and still a sufficient 
spectral resolution to distinguish absorption/emission features. 
The synthesized beam was set to 20", which is slightly larger 
than the best resolution achievable (~16"). The weighting pa¬ 
rameter was set to robust - 0.5, which is a combination of nat¬ 
ural and uniform weighting. These methods result in an rms of 
~14K (9mJy beam ') for emission free channels and up to a 
factor of 2 or 3 more in channels with strong emission, due to 
systematic errors such as side lobes. The dynamical range of our 
data set is ~ 100-200, depending on the region. In the follow¬ 
ing, the HI absorption and small scale structure is based on the 
THOR data, whereas the large scale emission is based on the 
VGPS data. 

The 21cm continuum data in this paper are taken from the Hi 
data cube for high and low velocities (-80 to -50 and 135 to 
155 km s') that are not affected by Hi emission or absorption. 
Therefore the data reduction and imaging for the 21cm contin¬ 
uum data is the same as for the H i data and we can avoid sys¬ 
tematic errors due to different data reduction methods. 


2.4. H 2 column density 


As we will compare the Hi and H 2 column density in Sect. 5.4 
we need an estimate of the H 2 content. We use dust observa 
tions from the Herschel Space Observatory to assess the H 2 col¬ 
umn density. These observations are based on the HiGal survey 
( Molinari et al.|2010]l. The H 2 col umn density map is taken from 
Fig. 9 of Nguyen-Lu’o’ng et al.| ( 2013|l, wh ich was derived by 
SED fitting methods described by Hill et al. ( 201 l| l. As the dust 
observations have no velocity resolution, we see all the dust and 
thus gas along the line of sight. |Carlhoff et aL] ( |2013| l showed that 
the Herschel dust data are similar within the uncertainties to the 
CO data at the velocity range of W43 (vlsr = 60 - 120 km s“'). 
Hence, the Herschel dust observations are dominated by the 
emission from W43 and the contributions of other regions along 
the line of sight can be neglected. We refrain from using CO 




Fig. 2. Sketch of the arrangement of the H i cloud with spin temperature 
Ts and optical depth t and the continuum source in the background. The 
on and off positions are marked. 


data to estimate the molecular content of W43, as CO does not 
trace all molecular hydrogen (e.g. |Smith et al.||2014] ) and CO 
becomes optically thick for the dense interior of W43 (jCarlhoff 
|etal.|2013^ . 


3. HI radiative transfer 


In this section, we explain the methodology we used to deter¬ 
mine the spin temperature, optical depth and column density of 
the neutral hydrogen. Even though the Hi 21cm line is a well 
known probe of these quantities, the arrangement of different 
Hi components with different temperatures along the line of 
sight can complicate these measurements. Eirst we will explain 
the classical methods to determine the column density via 
optically thin emission and Hi continuum absorption (HICA). 
Subsequently, we will outline the limitations of these methods 
and describe our approach to correct for the optical depth and 
the weak continuum emission. Although the basics of this 


approach are discussed in the literature (e.g., Wilson et al. 
|2^ , we modified the classical approach to account for the 
continuum emission and optical depth. Therefore, we outline 
the calculations here in more detail. 

All following values are frequency dependent. To keep the 
equations and description simple, we drop the frequency 
dependencies in the equations. We emphasise however, that 
the following emission/absorption mechanisms only work for 
identical frequencies and therefore identical velocities. 


3.1. Column density 


The classical approach to determine the Hi column density is 
given by [Wilson et al. ([2010 1 as: 


= 1.8224 X 10 '*' 

cm ^ 



( 1 ) 


where Vh, Ts and t are the column density, the spin temperature 
and the velocity dependent opacity, respectively. The spin tem¬ 
perature Ts describes the relative population of the spin states of 
the hydrogen atom ( [Wilson et al.|2010) l. As Ts is the equivalent 
of the excitation temperature for molecules, it is only equal to the 
kinetic temperature in Local Thermodynamic Equilibrium (i.e. 
when there are enough collisions to thermalise the gas). We as¬ 
sume in the following that the spin temperature T s does not vary 
significantly within one velocity channel. In most cases, neither 
Ts nor the optical depth are known. The simplest assumption is 
that the Hi emission is optically thin (see Sect. |3.2[ l. 
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Fig. 3. Optical depth simplification. The simplified optical depth is cal¬ 
culated by neglecting the spin temperature of the cloud. The solid and 
dashed lines represent a spin temperature of 50 K with optical depths of 
1 and 2, respectively. The dotted and dash-dotted lines represent a spin 
temperature of 100 K with optical depths of 1 and 2, respectively. 



r 


Fig. 4. The absolute uncertainty of the optical depth as a function of the 
optical depth itself for three different continuum background sources. 
The dashed, dotted and solid line represent continuum sources with a 
brightness temperature of Tcont = 500K, lOOOK and 2000K, respec¬ 
tively. The crosses on each line show the lower limit of the optical depth 
that we can observe for the given continuum brightness temperature. For 
larger optical depths, the absorption spectra saturates (see Sec. |4.1| for 
more details). 


3.2. Optically thin H i emission 

Under the optically thin assumption, without background con¬ 
tinuum emission, the expression for the brightness temperature 
Tb simplifies to the following; 


= Ts (1 - exp(-T)) ss Ts T. 


( 2 ) 


This simplification provides a linear relation between the col¬ 
umn density (equation[2l and the brightness temperature. Hence, 
we can estimate the column density directly from the measured 
brightness temperature T^. This method is used in several stud¬ 
ies and is well described in the literature (e.g., Lee et al.]|2012 


|Motte et al.|2014HWilson et al.|2010] l. Below, we will discuss its 
limitations and describe a procedure to overcome them. 
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3.3. Hi optical depth 


The Hi continuum absorption method is the classical method 
to determine properties of the cold neutral medium ( Strasser ^ 
|Taylor|[2004[ [Heiles & Troland]|2003a|bt [Strasser et al. 2007) 1. 
This method is based on strong continuum sources (Tb > 
300 K), such as Galactic Hii regions. Active Galactic Nuclei 
(AGNs) or extragalactic jets. As the brightness temperature of 
these continuum sources is larger than the typical spin temper¬ 
ature of the Hi clouds (7$ ~ 100K), we observe the Hi cloud 
in absorption. In addition, the absorption spectrum is dominated 
by the cold neutral medium, as the H i absorption coefficient is 
proportional to 

The classical observing strategy is the ‘on-off’ observation (see 
Fig .|^. The ‘on-source’ points toward the continuum source and 
reveals the absorption spectrum (Ton), whereas the ‘off-source’ 
points slightly offset from the continuum source and reveals the 
emission spectrum (Toff). During the data reduction, the contin¬ 
uum is not subtracted from the H i spectrum. Hence, we can mea¬ 
sure Tcont for channels that are not affected by the Hi line. The 
general assumption is that the on-source and off-source spectra 
originate from the same cloud with the same properties. There¬ 
fore it is important to have these two positions as close together 
as possible. 

The measured on-source and off-source brightness temperatures 
are: 


Ton = Ts (1 - exp(-T)) -H Toont exp(-T), 
Toff = Ts (1 - exp(-T)). 

Hence, the optical depth is: 



The spin temperature can then be determined from: 


(3) 


(4) 


Ts = 


Toff 


1 


(5) 


The advantage of the HICA method is the direct measurement of 
the optical depth and the spin temperature. The challenge of the 
HICA method is the need for strong continuum sources. Since 
most of the strong continuum sources are point sources, it is not 
possible to map the entire Milky Way, but instead the result is 
an incomplete grid of measurements. This way, it is not possible 
to map the intrinsic structure of individual clouds. Furthermore, 
the spin temperature does not describe the actual temperature 
of the cloud but is an mean of the cold and warm component, 
weighted by their column densities ( Strasser & Taylor)2004) 1. 
Consequently, the derived spin temperature is an upper limit for 
the cold component. 

As some Galactic continuum sources, such as W43-Main, are 
extended, it is difficult to determine a proper off position. Nev¬ 
ertheless, it is still possible to determine the optical depth. As 
we use the VLA Ch-D array observations including the contin¬ 
uum without the GBT data, we filter out most smooth large scale 
structure. The VLA C-hD array data show Hi emission of less 
than 30K. Therefore, we can neglect the emission of the H i cloud 
in equation and set Tog - (1 - exp(-T)) = 0 and calculate 

the optical depth without any off position as: 


simplified ~ 


(6) 
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Fig. 5. Ratio of observed and expected column density for a model cloud with spin temperatures of Tg = 50K and 100 K for the left and right 
panel respectively. The solid lines show the ratio for the assumption of optically thin emission. The dashed and dotted lines show the ratio for the 
assumption of optically thin emission and a diffuse continuum source with brightness 10 K and 30 K, respectively. 


Even without the effect of filtering out the H i emission due 
to the interferometer, this simplification holds for strong contin¬ 
uum sources (Tcont > 500K). Figure shows the relation be¬ 
tween the simplified and actual optical depth as a function of the 
continuum brightness temperature, not taking into account the 
filtering of the interferometer. It shows clearly, that even for high 
spin temperatures (Ts - 100 K) and high optical depth (r = 2), 
we miss less than 10% of the optical depth for a strong contin¬ 
uum source (Tcont ~ 2500 K). Therefore, we are able to measure 
the optical depth for the Hit region in the center of W43, even 
though we cannot determine a proper off position. 

We also consider the effect of optically thick clouds. For these 
optically thick clouds, the absorption spectra approaches zero. 
Due to the rms of the spectra, it is also possible that the absorp¬ 
tion spectra become negative, which is physically not meaning¬ 
ful. We set the optical depth to a lower limit for all absorption 
values that are close to zero. If the absorption spectrum Ton is 
smaller than 5 times the rms, we set the optical depth to: 


.„/5-tr(To„)\ 

^lower-limit —1^1 i. {/) 

\ cont / 

We will show that the optical depth saturates for the main 
velocity range of W43-Main, which has consequences for our 
interpretation and conclusions. 

The uncertainty of the optical depth (equation depends on 
Ton, Toff and Tcont- To estimate the uncertainty, we assume an 
uncertainty of 20K for all three quantities. The uncertainty of 
the optical depth is shown in Fig. as a function of the optical 
depth itself for three different continuum brightness values. It 
shows that the uncertainty increases significantly for increasing 
optical depths. However, we are not able to measure such high 
optical depths as the absorption spectra saturates for optically 
thick clouds and we report lower limits (see equation]^. These 
lower limits are marked with crosses in Fig. Up to the lower 
limit of the optical depth the corresponding uncertainty is ~0.3 
for the three different continuum brightness values. For strong 
continuum sources, such as W43-Main, the uncertainty of the 
optical depth is ~10% up to the lower limit of ~3. 


3.4. Column density corrections 

In addition to distinct, small continuum sources, such as Hit re¬ 
gions or extragalactic jets, we also find weak diffuse continuum 
emission in the Galactic plane. This component has a strength 
between 10 and 50 K. Therefore, it is not strong enough to in¬ 
duce absorption features (HICA), but nevertheless it can influ¬ 
ence the HI emission. To overcome this problem, the classical 
approach is to subtract this weak continuum emission during the 
data reduction. Hence, we subtract Tcont in equation and we 
get: 

Tb = Ts(l - eXp(-T)) -H Tcont eXp(-T) - Tcont 

= (Ts-Tcont) (l-exp(-T)). (8) 


This shows that even if we subtract the continuum from our data, 
it can still influence the measured brightness temperature. If we 
neglect the weak continuum emission, our measured H i emis¬ 
sion will be suppressed and thus the calculated column density 
will be underestimated. Therefore it is important to investigate 
the influence of the weak diffuse continuum emission. 

In the following, we will assess the assumption of optically thin 
emission and the influence of weak diffuse continuum emission 
on the determination of the measured column density for appro¬ 
priate model clouds. We assume two model clouds which have a 
spin temperature of 50 K and lOOK, respectively. We use equa¬ 
tion [T] to calculate the expected column density Nhi for differ¬ 
ent optical depths in one velocity channel {dv - 1.24kms *). 
Furthermore, we assume the cloud to be optically thin and use 
equation 1^ to calculate the brightness temperature of the Hi 
emission. Using this result, we can calculate the column den¬ 
sity of the cloud, but this time with the simplification of opti¬ 
cally thin emission. Hence, we call it the observed column den¬ 
sity N-niiobserved). The solid lines in Fig. [^show the ratio of 
the expected and observed column density for a range of optical 
depths. Obviously, for small optical depths (t < 0.1) our as¬ 
sumption of optically thin emission is sufficient and we observe 
more than 90% of the expected column density. But for larger 
optical depths (t > 0.1), we miss a significant fraction of the 
column density (>40%). 

In addition, we add weak continuum emission in the background, 
which changes the brightness temperature according to equation 
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[^and therefore suppresses the Hi emission. Nevertheless, if we 
still assume optically thin emission and do not consider the weak 
continuum emission, we can calculate the observed column den¬ 
sity. The dashed and dotted lines in Fig.|^show this case for con¬ 
tinuum emission of lOK and 30 K, respectively. Even for small 
optical depths (t < 0.1) we miss a significant fraction of the col¬ 
umn density, which depends on the ratio of the spin temperature 
and the continuum emission brightness temperature. In the worst 
case {Ts = 50K and Tcont = 30 K), we observe only 40% of the 
expected column density, even for small optical depths. 

To summarize, we measure the brightness temperature of the H i 
emission Tb, as well as the brightness temperature for the contin¬ 
uum Tcont and combine these information with the optical depth 
T, which we measure toward strong continuum sources. This al¬ 
lows us to calculate the corrected column density by combining 
equation and 

Nh = 1.8224 X 10'^ f_ (9) 

where the column density Nh is given in units cm“^, the bright¬ 
ness temperature Tb and the continuum brightness temperature 
Tcont are measured in K and the velocity dv is given in km s *. 


3.5. Continuum correction for strong point sources 

For the H i continuum absorption (HICA) toward strong contin¬ 
uum sources, we also have to consider weak diffuse continuum 
emission, which contributes to the on and the off position. There¬ 
fore we have to modify the equation for the on and off position 
(equation by adding another term for the weak diffuse con¬ 
tinuum emission, which we call Tcont,dif- This does not change 
our result for the optical depth (equation |^, as the optical depth 
depends solely on the difference between Ton and Toff. However, 
for the calculated spin temperature, we have to modify equation 
|5]to; 

rj. _ Toff “ Tcont,dif ^ /irw 


The effect of the diffuse continuum correction will be discussed 
in Sect. 15.21 


4. Results 


4. 1. Opticai depths determined using compact continuum 
sources 


As described in Sect. 3.3 we can use strong continuum sources 
in the background as light-houses, that shine through foreground 
HI clouds and create absorption spectra. Below, we characterize 
the continuum sources. 

To determine the optical depth of the Hi accurately, we need 
continuum sources in the background which are brighter than 
the typical Hi spin temperature (see Sect. |3.3| l. We extract 
all continuum sources in our field that have a brightness 
temperature Tcont > 400 K and a point-like structure, which 
yields in 8 point sources. The analysis of extended sources will 
follow in Sect. I4.4I We use a two dimensional Gaussian to ht 
the position and size of the continuum source. To measure the 
continuum brightness temperature, we average over high and 
low velocity channels in our Hi data cubes (-80 to -50 and 135 
to 155 km s '), that are not affected by the Hi line. The results 
are shown in Table [T] 

To determine the off spectrum (Toff), we average the emission 
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Fig. 6. H I emission and absorption spectrum of the extragalactic con¬ 
tinuum point source G31.388-0.384. The emission spectrum is shown 
in the first panel and is measured in an annulus around the point source 
with inner and outer radii of 60"and 120", respectively (corresponds 
to 3 and 6 times the restoring beam). The second panel presents the 
absorption spectrum toward the point source. The color shaded areas 
represent the approximate velocities of the Milky Way spiral arms jVaT| 
|lee|2008| ) in blue, red, green and yellow the Scutum-Centaurus, Cygnus, 
Sagittarius and Perseus arm, respectively, and the dashed black line in 
the first panel indicates the tangential velocity. In the first two panels, 
the grey shaded area indicates the 5cr noise level. The third panel shows 
the optical depth computed using equation]^ and the grey shaded area 
indicates the saturated optical depth limit, computed using equation]^ 
In the fourth panel the spin temperature is presented, which is computed 
using equation [T0| 


spectrum around the continuum source within an annulus 
with inner and outer radius of 60" and 120", respectively. 
The upper two panels of Fig. show a typical emission and 
absorption spectrum. We use equation]^ to calculate the optical 
depth for each channel. To avoid unrealistic optical depths, we 
only calculate the optical depth for those channels where the 
emission/absorption is 5 times larger than the corresponding 
noise. The grey shaded areas in the upper three panels in Fig. 
and|7]show the 5 sigma level. The optical depth is set to zero for 
channels, where the emission/absorption is below 5 times the 
corresponding noise. 

As explained in Sec. |3.3| we have to consider the effect of 
optically thick clouds for which we can determine only lower 
limits for the optical depth. The grey shaded areas in the third 
panels of Fig. and show these lower limits of the optical 
depth. 
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Table 1. Continuum point sources that were used for H i absorption studies. 


Name 

RA.(2000) 

(h : m : s) 

Dec.(J2000) 

MAJ X MIN 

(" X ") 

Tcont 

K 

^cont,dif 

K 

far Scutum 

arm 

spectral 

index 

location 

G31.388-0.384 

18:49:59.1 

-1.32.55.9 

21.3x20.3 

2055+24 

20.2 

yes 

-0.82 

ext. gal 

G31.411-h0.307 

18:47:34.1 

-1.12.45.5 

24.0 X 23.0 

914+26 

19.4 

no 

0.63 

gal. 

G30.234-0.138 

18:47:00.4 

-2.27.52.9 

21.7x20.7 

853+23 

42.1 

no 

0.00 

gal. 

G30.534-h0.021 

18:46:59.3 

-2.07.27.9 

28.0x21.5'* 

653+21 

35.9 

no 

0.21 

gal. 

G31.242-0. no 

18:48:44.7 

-1.33.13.7 

25.1 X 21.2 

573+22 

30.1 

no 

0.58 

gal. 

G29.090-h0.512 

18:42:35.8 

-3.11.04.7 

20.0 X 10.0 

571+44 

20.5 

yes 

-1.03 

ext. gal 

G30.744-H 1.008 

18:43:51.4 

-1.29.14.8 

23.1 X 20.1 

472+36 

16.9 

yes 

-0.90 

ext. gal 

G30.699-0.630 

18:49:36.3 

-2.16.27.2 

22.3 X 19.6 

403+18 

23.4 

yes 

-1.11 

ext. gal 


Notes. Continuum point sources in the observed field with brightness temperatures Tcont > 400 K. The names and spectral indexes are based on 
work presented in Johnston et al. (in prep) and are preliminary results. Coordinates, major and minor axes of the 2D fit, maximum brightness 
temperature T^oni of the continuum point source and the diffuse weak continuum temperature around the point source are given. The column ‘far 
Scutum arm’ indicates if absorption features at vlsr ~ -40kms“* are present in the Hi spectra. The spectral index and the absorption of the far 
Scutum arm are used to determine the location of the point sources, which is given in the last column. 

This source consists of two blended sources, which explains the large eccentricity. 



-50 0 50 100 150 

Vlsr [kms^^] 



Tspin [K] 


Fig. 8. Distribution of the spin temperature for all absorption features 
(655) toward point-like continuum sources. The median of the spin tem¬ 
perature is 97.5 K. 


cally thick regime. Nevertheless this overestimation is small, as 
for optical depths of t = 2 the measured spin temperature is 
Tspin ~ 1-16 ■ Toff (see equation!^. Therefore we overestimate 
Tspin by at most ~ 15%. We will mscuss these results in detail in 
Sect. IQ 


Fig. 7. Same layout as Fig. m but for the Galactic continuum source 
G31.411-1-0.307 


4.2. HI spin temperature toward compact continuum sources 

For each channel that allows for an optical depth measurement, 
we compute the spin temperature using equatio n [TO] The spin 
temperature is shown in the fourth panel of Fig. [6[and [7] This 
method reveals absorption features in 655 channels for all eight 
continuum sources. The median of the spin temperature is 97.5 K 
and the distribution of all absorption features is shown in Fig. 

For large optical depths, the spin temperature approaches the 
brightness temperature of the off-position. However, for large 
optical depths (t > 2) we can only report lower limits for 
r. Hence, we overestimate the spin temperature in the opti- 


4.3. Location of continuum sources 


To characterize and understand the Hi absorption spectrum 
toward the continuum point sources, it is important to know the 
location of the continuum source. We can distinguish between 
extragalactic and Galactic point sources. The Galactic sources 
are mostly H ii regions, whereas the extragalactic point sources 
are radio lobes from extragalactic jets or AGN. To distinguish 
between them, we have two criteria: the spectral index and H i 
absorption of the far Scutum-Centaurus arm. As we have to 
consider e.g. primary beam effects or different spatial filtering 
of the interferometer for different frequencies, it is very difficult 
to determine accurate spectral indexes (e.g. |Rau et al.||20T4l 
|Bhatnagar et al. [2011 |2013p . For this analysis, the spectral 
indices were calculated using flux measurements in the two 
most-separated spectral windows at 1.05 and 1.95 GHz and are 
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Fig. 9. 1.4 GHz continuum emission of W43-Main. The yellow star in¬ 
dicates the position of the OB-cluster ( [Lester et al.|1985| > and the white 
crosses marks the position of UCH it regions observed by the CORNISH 
survey ( jHoare et al. [2012] [Purcell et al.|20l3t 


brightness temperature for this region is Tcont ~ 2200 K and 
therefore we are able to calculate the optical depth, even 
though we cannot determine a proper off position (see Sect. 
[3.3| for details). The Hi absorption spectra shows a cut-off at 
visr ~ lOOkms^' (see Fig. 10 1 , which marks the velocity of 
the continuum source. We measure radio recombination line 
emission at the same velocity (Johnston et al. in prep, [Anderson] 
[et al.[|^ll| l. Therefore, the continuum source W43-Main is 
situated at visr ~ 100km s“\ which is close to the tangent point 
velocity of the Scutum-Centaurus arm (jNguyen Luong et al.[ 

[20TTI 1. 

Furthermore, the observed field contains two other extended 
continuum sources: W43-South and the supernova remnant 
SNR G029.7-00.2 with continuum brightness temperatures of at 
most Tcont ~ 1840 K and Tcont ~ 850 K, respectively. 


based on work presented in Johnston et al. (in prep.). 

Four sources in our sample have negative spectral indices (Table 
[TJ, which is typical of synchrotron radiation from extragalactic 
jets or AGNs. The remaining four sources show flat or posi¬ 
tive spectral indices, which indicates free-free emission from 
potentially optically thick Hii regions. These four free-free 
emission sources also show no Hi absorption feature for the 
far Scutum-Centaurus arm. Therefore, these sources reside in 
the Milky Way. Furthermore, we can study the Hi absorption 
spectra and molecular emission spectra of the Galactic point 
sources to estimate their distance and overcome the near-far 
distance ambiguity. This is done by Anderson et al. ( 2014| l 
and they find that e.g. the source G31.411+0.307 has a near 
distance of 6.6 kpc. For this source we also see a sham cut off 
in the absorption spectra at visr ~ 100 km s“' (see FigjTkwhkh 
corresponds to the molecular line velocity reported in (Anders on] 
[etar] ( |20T4) . 

Figure shows an example of an extragalactic continuum 
source. The characteristic absorption of the far Scutum- 
Centaurus arm at visr a; -40kms“' is clearly observable. 
Furthermore, the absorption and emission spectra show a similar 
cut off for high velocities, when approaching the tangential 
point of the near Scutum-Centaurus arm. On the other hand, 
the absorption spectrum of Fig. [^neither shows the absorption 
of the far Scutum-Centaurus arm, nor approaches zero at the 
velocity of the tangential point. This is typical for Galactic 
continuum sources. 


4.5. Optical depth of W43-Main 


In Sect. [3.3[ we described that we can determine the optical 
depth for strong continuum sources by using the absorption 
spectrum. As the brightness of the continuum source W43-Main 
is Tcont ~ 2200 K, the uncertainty for the optical depth measure¬ 
ment is ~ 10% (see Fig. [^. Figure [TOj shows the optical depth 
of W43-Main and the grey shaded area indicates the lower limit 
of our measurement with Tiower-iimit = 2.7. The optical depth 
peak at visi ~ 10 km s“' can be allocated to the near Sagittarius 
arm and therefore is not connected to the actual star forming 
region W43. In contrast to that, it is difficult to allocate the 
distinct absorption features between 50 and 80 km s ' and it 
is not clear whether they are spatially connected to W43. The 
prominent absorption feature of W43 is situated between 80 and 
100kms '. In this region, our measurement is saturated and 
therefore we can only report lower limits for the optical depth. 
As we resolve the strong continuum source W43-Main spatially, 
we can determine the optical depth along different lines of 
sight and thus investigate the spatial distribution of absorption 
features as done in detail by jLiszt et al. ( 1993| l. We refrain from 
such a study as we are mostly interested in the velocity range of 
W43 (visr = 80 - llOkms^') at which the absorption features 
are mostly saturated preventing a detailed study of the spatial 
distribution. We instead measured the optical depth toward the 
strongest continuum peak to estimate the maximum optical 
depth possible (see Fig. [TOli, which is nevertheless a lower limit. 


4.4. Extended continuum sources - W43-Main 

Besides point-like continuum sources, our observed field 
contains three strong {T > 800 K) extended continuum sources. 
The most prominent source is the well known H ii region around 
W43-Main ( Lester et al.||1985 1 , which is shown in Fig. This 
region has an angular extent of ~ 300", which corresponds to 
~ 8 pc at a distance of 5.5 kpc. As our resolution is 20", we are 
able to resolve the internal structure well. It is known that the 
nebula is ionized by an OB-cluster and the observed continuum 
signal is the result of free-free emission. Figure shows the 
continuum emission of W43-Main and the yellow star marks 
the position of the OB-cluster (J2000, 18h47m36s, -1°56'33", 
[Lester et al.| 1985| l. Furthermore, several UCHii regions can be 
found in the CORNISH survey, which are marked with white 
crosses ( [Hoare et al.[[2012] [Purcell et al.[[2013| l. The maximum 


4.6. Hi column density with optically thin assumption 

In this section we assess the Hi column density and the Hi 
mass of the GMC associated with W43. These two quantities 
depend strongly on the chosen velocity range. [Nguyen Luong[ 
[et al.[ ( |201 Ij l defined for the ‘main’ and ‘complete’ velocity 
range of W43 values of vx^^fmain) = 80 - llOkms^* and 
visAcomplete) = 60 - 120kms“*. In the following, we will 
use the ‘complete’ velocity range, but we stress that a different 
velocity range can change the column density and mass signih- 
cantly. 

As described in Sect. [3.2[ we can determine the column density 
by assuming optically thin Hi emission. For this method, the 
column density is proportional to the observed Hi brightness 
temperature. The column density map assuming optically thin 


emission is shown in the top-right panel of Fig. 11 To estimate 
the mass, we integrate this column density over the main region 
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Fig. 10. Optical depth of W43-Main, calculated using equation]^ The 
grey shaded area indicates the maximum measurable optical depth of 
Tiower-iimit = 2.7 Calculated using equation]?] 


of W43 (I - 29.0 - 31.5°, \b\ < 1°). We mask all regions, where 
the emission spectrum has negative values, and exclude them. 
Given a distance of 5.5 kpc (]Zhang et al.]]2014)l, our analysis 
finds an Hi mass of M = 2.7x 10 “Mq (/ = 29.0-31.5°, \b\ < 1°, 
v’lsr = 60- 120 km s“'). 


Weak diffuse continuum emission can influence the H i emission 
spectrum and therefore the H i column density calculation needs 
to be modified according to equation]!] (Sect. ]T4) i. The top-right 
panel of Fig. 11 shows this effect clearly. The color represents 
the column density determined by the optically thin assumption 
and the black contours indicate the weak diffuse continuum 
emission. These two components show a clear anti-correlation. 
However, this anti-correlation is the result of the expected ob¬ 
servational effect in which H i emission is suppressed by weak 
continuum emission (see Sec. ]3.4| l. To overcome this problem 
we have to consider the optical depth, which we discuss in the 
next section. 


4.7. Hi column density with optical depth correction 


With the measurement of the optical depth and the weak contin¬ 
uum emission, we can correct the H i emission as explained in 
Sect. 13.41 This allows a more accurate determination of the col¬ 
umn density, which is shown in the middle panel of Fig. [m We 
chose the ‘complete’ velocity range of visr = 60 - 120km s“'. 
As we correct for the optical depth, we are able to observe a 
larger Hi column density. We also correct for the weak con¬ 
tinuum emission around W43-Main. This correction removes 
the anti-correlation between the continuum emission and the H i 
column density. Hence, we determine the H i mass to be M = 
6.1 X 10® Mo (I = 29.0 - 31.5°, \b\ < 1°, v^r = 60- 104kms-'). 


The optical depth spectrum of W43-Main (Fig. ]T0| | shows, that 
the absorption ends abruptly at lOOkms \ but the Hi emission 
as well as C O emission ( jNguyen Luong et al. 2011 Carlhoff 


et al. 


2013 I reveals features up to visr = llOkms V The rea¬ 


son for the abrupt drop in the absorption spectrum is not the 
absence of Hi, but the location along the line of sight of the 
continuum source at this velocity (see Sect. ]4.3] for further de¬ 
tails). As we do not see H i in absorption for velocities larger than 
100kms ', we are also not able to determine the corresponding 
optical depth. Therefore we cannot apply our corrections to the 


Hi column density for velocities larger than lOOkms^*. Hence, 
the velocity range for the previously mentioned H i mass is only 
Visr = 60 - 104kms“*. Nevertheless we can determine the Hi 
mass for the velocity range from visr = 104 - 120kms“* by 
using the optically thin assumption. This reveals a Hi mass of 
M - 0.5 X 10® Mq. Hence, the total Hi mass for W43 in the 
velocity range visr = 60 - 120kms“* is M - 6.6 x 10® Mq 
with the optical depth correction for the velocity range vis, = 
60 - 104 km s ' and optically thin assumption for the velocity 
range visr = 104 - 120kms“'. This is 2.4 times larger, than the 
H I mass determined with the optically thin assumption. 

The limitations and uncertainties of our determined H i column 
density and Hi mass with the optical depth corrections will be 
further discussed in Sect. 15.31 


4.8. Hi self absorption 

The Hi Self Absorption (HISA) method uses the diffuse broad 
Hi emission background of the Milky Way as illumination for 
colder foreground clouds (e.g. ]Gibson et al.]]2000 2005a|b O] 


]& Goldsmith|2003]]McClure-Griffiths et al.|2b06 l. Dark regions 

on maps and narrow absorption features in spectra reveal these 
HISA features. The terminology ‘self absorption’ can be mis¬ 
leading, as the H i emission and absorption can occur in the same 
cloud or at the same position but does not have to. The advan¬ 
tage of this method is that it is possible to map entire absorption 
clouds and study their intrinsic structure. On the other hand, the 
disadvantage is that a sufficient background emission with the 
same velocity as the absorbing cloud is necessary to detect HISA 
features. Therefore this method misses a large portion of the cold 
HI clouds. The differentiation between actual HISA features and 
the lack of H i emission can also be challenging. ]Gibson et al. 
( 2005bj l, however, present an efficient method to detect HISA 
features automatically. Another disadvantage is that the optical 
depth and spin temperature can only be measured together and 
further assumptions are needed to disentangle these two values. 
As described in the previous section, we correct for the opti¬ 
cal depth effects and weak continuum emission. This correction 
does not account for locally confined HISA features, as we as¬ 
sume a uniform optical depth for the entire W43 region. HISA 
features could have a higher and spatially varying optical depth 
that we cannot measure. Furthermore, the weak and diffuse con¬ 
tinuum emission around W43 makes the search for HISA fea¬ 
tures even more complicated. Hence, we refrain from searching 
and analyzing the possible HISA features around W43. There¬ 
fore we are likely missing some cold neutral hydrogen in our 
analysis. HISA studies will be conducted in other parts of the 
Milky Way THOR study with less diffuse continuum emission. 


5. Discussion 


5.1. Phases of the Neutral Atomic Hydrogen 

It is well known that the neutral atomic ISM has several phases 
that coexist side by side with very different properties (e.g.]Clark 
IWolfire et al.|]T995l [Heiles & Troland]|2003bl |Wolfire 


1965 


et al. 


20031. The main constituents are the cold neutral medium 


(CNM) and the warm neutral r nedium (WNM) with spin temper¬ 
atures on the order of <100K ( Strasser et al.|20()7 1 and ~ 10‘'K 
( JMurray et al. 2014] Roy et al. 2013a i, respectively. Further¬ 
more, their density differs by two order o f magnitude (CNM : 


Uh ~ 50cm WNM: Uh ~ 0.5cm e.g.,]Stahler et al.|2005 i. 


Due to the different spin temperature and density of the CNM 
and WNM, their corresponding optical depths are significantly 
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Fig. 11. The top-left panel shows the continuum emission at 21cm in Kelvin. The top-right and middle panels presents the Hi column density 
for optically thin assumption and optical depth corrections, respectively. The continuum and Hi emission data are based on VGPS, whereas the 
optical depth correction used in the middle panel is based on the THOR data. The bottom panel shows the H 2 column density from the HiGAL 
data ( [Nguyen-Lu’o’ng et al.|201^ . In all panels, the black and white/blue contours present the continuum emission at 21cm (black contours show 
levels of 10, 30 and 70 K, white/blue contours show levels of 200, 400, 600 and 800 K). The black ellipses in the middle and lower panel have an 
equivalent radius of r = 80 and 140pc. Several important objects are marked. 
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different with typical values of twnm ~ 10^-10“* (Murray 
et al.||2014| l and tcnm ^0.1 ( Strasser & Taylo^|2004| l. This is 


important for our interpretation. 

Due to the low optical depth of the WNM, we see in absorp¬ 
tion merely the CNM. Hence, the optical depth shown in Fig. 


10 is the optical depth of the CNM. As the absorption spectra 


are strongly dominated by the CNM, we are not able to measure 
the optical depth of the WNM and we can assume the WNM 
to be optically thin. In contrast to this, we see a combination of 
the CNM and WNM for the H i emission. For the correction of 
the column density (see Sec. 4.7 1 we use the optical depth in¬ 
formation from the absorption study to correct the H i emission 
data. As we do not distinguish between the two phases, we might 
combine two different quantities, namely the optical depth of the 
CNM with the emission of the CNM and the WNM. This could 
lead to an overestimation of the column density. In the follow¬ 
ing, we will assess this effect. 

If we assume a CNM cloud with varying Tspin(CNM)~20- 
80K and varying optical depth surrounded by a WNM with 
Tspm(WNM)~7000K (Murray et al. 2014 1 and optical depth 
t(WNM)~5x 10“^ we can calculate the column density of each 
component separately (CNM and WNM). Furthermore, we can 
calculate the brightness temperature that we would observe and 
apply our correction method described in Sec. |3.4| Finally, we 
can compare the actual column density of the CNM and WNM 
with the column density we would observe with our correction 
method and investigate the overestimation of the actual column 
density. The result is shown in Fig. which shows the ratio of 
the observed and actual column density as a function of the op¬ 
tical depth of the CNM t(CNM) for different CNM spin temper¬ 
atures Tspin(CNM). Fig. 12 shows that we measure the column 
density accurately for the optically thin case and we overestimate 
the column density for larger optical depths. However, even for 
the extreme case (Tspin(CNM)~20K, t(CNM)~ 3), we overesti¬ 
mate the column density by at most 1.35. This effect is smaller 
than the underestimation of the column density due to saturated 
optical depths and therefore we will consider only a single cold 
component in the following. Similar results in simple simula¬ 
tions for different combination of CNM and WNM fraction and 
a wide range of N// and t were found by Chengalur et al. ( |2013| l 
and |Roy et al.| ( |2013a|b| l. 


5.2. HI spin temperature measurements toward point 
sources 

The spin temperature is the best measure of the kinetic tempera¬ 
ture of the HI gas along the line of sight. Hence, it is a principal 
determinant of the physical processes that occur in the H i gas. 
As described in Sect. |3.5| we corrected the H i spin tempera¬ 
ture for diffuse weak continuum emission. If we neglect this 
correction we would systematically overestimate the Hi spin 
temperature and our sample would have a median Hi spin 
temperature of llOK. That means the median value would be 
~ 12 K higher compared to our corrected value of 97.5 K (see 
Fig.§. 

Similar studies of the spin temperature with similar methods 
can be found in the literature. For example [Strasser & Taylor] 
( |2004| l report a median spin temperature of 120 K. This is 
~ 20 K higher than our value of 97.5 K, but they do not consider 
the diffuse weak continuum emission and therefore probably 
overestimate the spin temperature. 

[Heiles & Troland] (|2003a|b|l developed an extensive method to 



Fig. 12. The ratio of the observed and actual column density as a func¬ 
tion of the CNM optical depth t(CNM) for different CNM spin temper¬ 
atures Tspi„(CNM). 


fit Gaussian components to the absorption and emission spectra. 
Spectral features shown in absorption and emission are assigned 
to the cold neutral medium (CNM), whereas emission-only 
features are assigned to the warm neutral medium (WNM). 
Using this method, they are able to distinguish these two phases 
and measure their properties such as as the spin temperature 
individually. However, they report that it is difficult to use their 
method for sources close to the Galactic plane {\b\ < 10 deg), as 
multiple components can overlap and the Galactic rotation can 
broaden their profiles. Figure]^ and [7] illustrate this problem for 
our region and for these spectra it is impossible to find unique 
Gaussian components. As |Heiles & Troland] ( |2003a|b| l are able 
to fit individual components, they find on average lower spin 
temperatures for the CNM in the range of ~ 40 - 70 K. Given 
that we can not distinguish between the CNM and the WNM, 
the measured spin temperature is a column density weighted 
mean between these two components ( [Strasser & Taylor|2004) l. 
This could explain the increased spin temperature reported in 
Sec. 4.2 Another explanation for the increased spin temperature 
could be the strong radiation field of W43. 

Further complications in our sample are the Galactic continuum 
sources. As these only trace the optical depth up to their location 
within the Milky Way, we miss the optical depth of the hydrogen 
behind the continuum source. On the other hand, the off-position 
measurement traces all the hydrogen along the line of sight. As 
the kinematic distance is uncertain due to near/far ambiguities, 
we cannot use the velocity to distinguish the distance of the 
emitting and absorbing hydrogen. Hence we underestimate the 
optical depth for the Galactic continuum sources and therefore 
overestimate the spin temperature. However, we do not see a 
significant difference in the mean spin temperature for Galactic 
and extragalactic continuum sources and we can neglect this 
effect. 


5.3. Coiumn density and mass estimate 

In Sect. [4.7[ we explained that we correct the H i column den¬ 
sity for optical depth effects as well as for the diffuse continuum 
emission. This leads to a more accurate estimate of the column 
density and a more accurate mass estimate. Nevertheless these 
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corrections have limitations which we discuss in the following 
section. 

As we need a strong continuum source in the background to de¬ 
termine the optical depth, we can measure the optical depth only 
toward certain locations. For our Hi column density correction, 
we used the strong continuum source W43-Main. Furthermore 
we assumed that the optical depth is the same for the entire 
cloud. This assumption might not hold, especially for the outer 
parts of the cloud. To investigate this effect, we compared the 
optical depth measurements for other sources at other positions. 
However, in the range of visr = 80 - 110 km s ', we find that the 
optical depth is mostly saturated and determined by the corre¬ 
sponding lower limit. Two examples are given in the third panel 
of Fig. and which reveal a lower limit of riower-iimit = 2.9 
and Tiower-iimit = 1 -9, respectively. Other continuum sources that 
have a larger separation from the Galactic mid plane show sim¬ 
ilar results. For example the continuum source G30.699-0.630 
has a Galactic latitude of b ^ 0.6°and still shows a saturated op¬ 
tical depth for the velocity range of W43 with a lower limit of 
Tower-limit = 1.5. If we use this continuum source to correct for 
optical depth effects and the weak diffuse continuum emission, 
we determine a mass of M = 4.8 x 10^ Mq for the whole cloud 
{I = 29.0 - 31.5°, \b\ < r, Vkr = 60 - 120km s^'). Hence, by 
using W43-Main to correct the optical depth for the entire cloud, 
we do not overestimate the mass in the outer parts significantly, 
but at most by a factor of 1.4. We will also use this mass estimate 
as a lower limit for the mass of W43. 

The saturation of the optical depth is the second limitation we 
have to consider, especially for the inner part of the cloud. In 
this region we underestimate the opacity and, therefore, the col¬ 
umn density which would lead to a further underestimate of the 
mass. Estimating this effect is difficult and therefore we report 
only lower limits for the column density and the mass in the in¬ 
ner part of W43. 

The third limitation is the distance of the continuum source that 
we use to determine the optical depth. As explained in Sect. |4.3[ 
the continuum source W43-Main is located close to the tangen¬ 
tial point of the Scutum-Centaurus arm at a distance of 5.5 kpc. 
Hence, we only see Hi absorption as far as this distance, and 
miss all the H i that is located behind the continuum source but 
still within the cloud. If we assume that the continuum source is 
at the center of the cloud, we underestimate the column density 
by a factor of two. Another approach to overcome this limita¬ 
tion is to look at other more distant continuum sources. The con¬ 
tinuum source G31.388-0.384 is extragalactic and has a com¬ 
parable brightness temperature to W43-Main. Hence, it is an 
ideal candidate for this test. We used the optical depth shown 
in the third panel of Fig. to correct the H i column density 
and measured the corrected mass for the same area and veloc¬ 
ity range, as in Sect. |4.7| The absopition spectrum of W43-Main 
and G31.388-0.384 are similar, except for the velocity range 
Hsr = 100 - 120 km s“\ which is missing in the W43-Main 
spectrum. We can also use the optical depth measurements of 
G31.388-0.384 to correct the Hi emission. However, as these 
absorption spectra are similar, the corrected masses are the same 
within the uncertainties (M = 6.9 x 10^ Mq for the correction 
with the optical depth of G31.388-0.384). 

Summing up, we might overestimate the mass in the outer part 
of W43, but underestimate the mass in the inner part of W43. 
Due to the position of the continuum source, we might underes¬ 
timate the mass by a factor of two. As explained in Sect. |4.8| we 
also miss some H i due to self-absorption. Hence we report the 
H I mass of W43 to be a lower limit of M - 6 . 6 - 1.8 x 10^ Mq. 
Several Hi mass estimates are given in the literature (|Nguyen| 
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Fig. 13. The column density of H i and H 2 measured in elliptical annuli 
around W43 MMl. The x-axis presents the equivalent radius of these 
elliptical annuli. The black and green lines represent the fitted curves 
with the corresponding Icr uncertainties shown as a grey shaded area. 
The larger symbols (diamonds for H 2 , squares for Hi) present the av¬ 
eraged value of the elliptical annuli and their corresponding variations 
within the annuli. As the optical depth spectra saturates, we can only de¬ 
termine lower limits for the H i column density. The blue dots and line 
show the H 2 /HI column density ratio of the data and the fitted curves, 
respectively. 


|Luong et al.|201 l[|Motte et al.|2014| l. All of them are calculated 
with the assumption of optically thin emission. 

Motte et al.|(|20141l measured the Hi mass in the inner part of 
W43 {I = 29°.6 ~ 3r.4 and b = -0°.5 ~ 0°.3) and for the ve¬ 
locity range (visr = 60- 120kms“'). The assumption of optically 
thin emission reveals an H i mass of M - 0.9 x 10® Mq. If we 
use our corrected H i column density to measure the mass in this 
area, we determine a H i mass of M - 2.2 x 10® Mq (this value is 
smaller than the value given in Sect. |4.7| as we only consider the 
inner part of W43). As outlined above, we claim that this mass 
estimate is a lower limit, and therefore the mass determined with 
the optically thin emission is at least a factor of 2.4 too small. 
This has implications for theoretical models, which we will dis¬ 
cuss in Sect. 15 : 6 ] 


5.4. Spatial distribution of hydrogen 

Beside the total H i mass, knowing the spatial distribution of the 
Hi is crucial to understand the formation of W43. As we have 
corrected the column density map for optical depth effects and 
the weak diffuse continuum emission, we can use this corrected 
column density map (see Fig. [n) to investigate the Hi spatial 
distribution in more detail, especially in the center. 

Similar work was done by |Motte et al.| ( |2()14| l. They measured 
the H i column density in rectangular annuli around the center of 
W43 (/ = 30.5°, b = 0°) with an aspect ratio of 3:2 and find an 
increasing H i column density inwards from Ahi ~ 40 Mq pc“^ at 
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Fig. 14. The density of H i and H 2 as a function of the distance toward 
the center of W43 (equivalent radius). The diamonds and the dashed 
line represent the molecular hydrogen, whereas squares and the solid 
line represent the atomic hydrogen. The corresponding 3cr uncertainties 
from the fit introduced in equation[T^are shown as grey shaded areas. 


a distance of 170pc to A^hi ~ 80Mopc“^ at a distance of 50pc 
from the center (velocity range visr = 60 - 120kms“'). Since 
they assume optically thin H i emission, they underestimate the 
HI column density, especially in the central region (as discussed 
in Sect. |5.3| l. 

To compensate for the approximate elliptical structure of W43, 
we choose elliptical annuli with an aspect ratio of 3:2 for ma¬ 
jor and minor axis, which fits the large scale structure of W43 
well (see Fig. As we focus on W43-Main, we choose the 
most massive submillimeter peak W43-MM1 {I - 30.8175° and 
b = -0.0571°, |Motte et al.|2003] l as the center for the ellipses. 
The width of each annulus is 10 pc for the major axis and 6.6 pc 
for the minor axis. For each elliptical annulus, we calculated the 
equivalent radius r = ^Jma jor ■ minor and assigned these values 
as the distance to the center shown in Fig. [1^ 

The black squares in Fig. represent the H i column density 
mean value and the corresponding standard deviation of each el¬ 
liptical annulus. We conhrm the result of |Motte et al.| ( |2014| l, 
that the Hi column density rises inwards. Our corrections allow 
us for the first time to study the central region of W43 (r < 50 pc) 
and, within the uncertainties, we report a flat column density dis¬ 
tribution in this inner region. However, we mention that this flat¬ 
ness could also be due to the underestimation of the column den¬ 
sity in the central part as the optical depth saturates and therefore 
the column density is a lower limit. 

Furthermore, the diamond symbols in Fig. [T^ present the col¬ 
umn density of the molecular hydrogen. The H 2 distribution is 
centrally concentrated and the column density rises steeply to¬ 
ward the center, which is different from the Hi profile. The 


large uncertainties in the center are due to the dumpiness of 
the molecular hydrogen as the most prominent clumps, such 
as W43-MM1, are located in the first two elliptical annuli. Be¬ 
side W43-Main, the second most prominent molecular clump is 
W43-South (see Fig. 111. However, in this analysis we focus on 
W43-Main and therefore, we choose W43-MM1 as the center for 
the ellipses. Furthermore, we masked W43-South for this anal¬ 
ysis, as it would introduce a large uncertainty at an equivalent 
radius of r ~ 50 pc. 

We use a non-linear least square method (’curveflt’ in the scipy 
package) to fit the mean values of the thin elliptical annuli with 
their corresponding uncertainties. For the H i column density we 
assume a quadratic radial distribution; 


Nniif) - a ■ r + b ■ r + c. 


( 11 ) 


where A^hi(^) describes the Hi column density and r describes 
the equivalent radius. The free parameters a, b and c have the 
fitted values of a = -3.9 + 0.7 x 10“^ Mopc^^, b - 0.089 + 
0.099 Mq pc“^ and c - 171 + 3Mopc“^. For the H 2 distribution 
we assume an exponential function to fit the data; 


Nn-,(r) = d ■ exp(-e ■ r) + f. 


( 12 ) 


where A^H 2 (^) describes the H 2 column density and r describes 
the equivalent radius. The free parameters d, e and / have the 
fitted values of d - 458 + 44Mopc“^, e - 0.022 + 0.004pc“' 

and/ = 97±15Mopc-2. 

The black and green lines in Fig. [T^ present the column den¬ 
sity fits and the grey shaded areas represent their correspond¬ 
ing uncertainties. Both, the Hi and H 2 distributions are well 
fit by the assumed functions. As the uncertainties for the in¬ 
nermost part of the H 2 distribution (r < 20 pc) are large, the 
fitted curve deviates from the data points and this area has to 
be treated cautiously. We also tested different functions to fit 
the data, such as functions with only two free parameters, or a 
Gaussian distribution for the H i distribution, but the results were 
similar within the uncertainties. Finally, we chose the aforemen¬ 
tioned functions as they could reproduce well the ratio of the H 2 
and Hi column density. This ratio is shown in blue in Fig. 13 
The circles represent the ratio of the data points, whereas the 
solid line represents the ratio of the fitted curves. The ratio stays 
fairly constant at ~ 1 for 140 pc > r> 60 pc. For this re¬ 

gion, the HI column density also rises to its maximum value of 
Ahi ~ 170Mopc“^. Further inwards (r < 60pc), the Hi column 
density stays constant at this maximum value, whereas the H 2 
column density rises sharply. Hence, the H 2 /H i ratio also rises 
sharply to values above three. Summing up, the column density 
measurements imply that we have a mixture of H i and H 2 in the 
outskirt of the cloud and a molecular dominated region in the 
center. In the following, we will investigate this structure for the 
particle density. 

As mentioned before, we And a flat column density for the H i 
distribution toward the center of W43-Main, but what does this 
imply for the actual density in the center? If we assume that 
W43 has an elliptical shape, we can decompose the cloud into 
different layers, similar to an onion. Furthermore, we assign the 
appropriate column density to each layer with the information 
given in Fig.[T^ As the column density is additive, the appropri¬ 
ate column density of each layer is the measured column density 
at the position of the considered layer minus the column density 
of all layers further outside. Hence, a flat H i column density dis¬ 
tribution toward the center (see Fig.[T3]l means that the layers in 
the center have no or a very small column density and therefore 
also a very small particle density. In the following section, we 
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will use the concept of a elliptical layered structure to determine 
the actual particle density and show that indeed the H i particle 
density drops toward the center of W43 within this model. 


5.5. Linking coiumn density to particie density 


While the column density neglects the third dimension, it does 
not necessarily reflect the actual particle density. Hence, we have 
to estimate the third dimension. Using the elliptical, layered 
structure for W43 previously described and the results presented 
in Fig.[T^ we estimate the particle density in each layer; 


^ (M-M+i) 

2 (n+i - n) ’ 


(13) 


where tii describes the particle density of layer i, A^i and A^i+i 
describes the column density of layer i and layer i+1, respec¬ 
tively, and Ti and r-^+i describes the equivalent radius of layer i 
and i-Hl. The factor of two accounts for the two layers of the el¬ 
liptical annuli in the front and the back of the cloud. The result 
for this calculation is shown in Fig. as squares and diamond 
data points for the H i and H 2 density, respectively. To increase 
the signal to noise ratio, we used larger elliptical annuli with a 
major axis of 20pc and a minor axis of 13.3 pc. 

In addition, we can use the fitted curves of the column density 
to estimate the density as well, by inserting equation [TT| and [T2] 
into equation[T^ Hence, the radial distribution of the H i and H 2 
densities are; 


mhi = -0.5 (a (n H- n+i) - b), 

«H 2 = T , -r (exp(-e • n) - exp(-e ■ n+i)), (14) 

2 (n+i - n) 


where hhi and hhi describes the particle density of Hi and H 2 , 
respectively, ri describes the equivalent radius of layer i and a, 
b, c and d are the free fit parameters introduced in equations 
In Fig. 14 the calculated density distributions are 


11 and 12 


shown as black lines with the corresponding uncertainties as a 
grey shaded areas. We note that the density distributions of Hi 
and H 2 are very distinct. While the H i distribution follows a sim¬ 
ple linear relation, the H 2 distribution shows an exponential in¬ 
crease toward the center. The particle density shows a mixture 
of the atomic and molecular hydrogen in the outskirt of W43 
(r > 100 pc) and a molecular dominated interior, similar to the 
measurements of the column density. However, the particle den¬ 
sity of the atomic hydrogen drops to almost zero toward the cen¬ 
ter which results in the observed constant column density. 
Another way to present our results is shown in Fig. This fig¬ 
ure shows the particle density of H i and H 2 introduced in equa¬ 
tion 14 as a function of the total column density (Nhi -i- Nhj). 


High column density regions, i.e. the central region of W43, are 
dominated by molecular hydrogen. On the other hand, we find 
an equivalent mixture of atomic and molecular hydrogen for low 
column density regions, which represent the envelope of W43. 

In addition, we study the correlation between the Hi and H 2 
density. To do so, we replace the equivalent radius in Fig. 
and plot the H 2 density as a function of the H i density. Figure 


16 shows the corresponding plot with the uncertainties for the H 2 


density. The statistical and systematical uncertainties for this plot 
are relatively large, but nevertheless we see an anti-correlation 
between the Hi and H 2 density. This can be explained with a 
simple model that H 2 forms out of H i in the innermost part of 
W43. However, we do not detect a sharp transition between H 2 
and H i predicted by Krumholz et al. (2008 2009) 1. In the follow¬ 
ing section, we will discuss possible implications. 
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Fig. 15. The particle density of Hi and H 2 as a function of the total 
column density (Nh, + Nhj). The diamonds and the dashed line repre¬ 
sents the molecular hydrogen, whereas squares the solid line represents 
the atomic hydrogen. The corresponding 3<t uncertainties from the fit 
introduced in equations[TT]and[^are shown as grey shaded areas. 



Fig. 16. The particle density of molecular hydrogen as a function of the 
atomic hydrogen. The black curve represents the best fitted curves from 
Fig. and the grey shaded area shows the 3cr uncertainty for the H 2 
fit. 

5.6. Coiumn density threshoid for atomic hydrogen 

A fundamental question of molecular cloud formation is how 
does molecular hydrogen form out of neutral atomic hydrogen 
and what are the corresponding conditions ( |Dobbs et akpOlTl l. 
It is well known that the density must be high enough to shield 
the molecular hydrogen from the interstellar radiation field, to 
avoid dissociation back to its atomic form (|Hollenbach & Tie-) 
liens) 1997) [McKee & Ostriker|20071 [Krumholz et al.|2008| 12009) 
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Table 2. Densities and model parameter. 


r 

[pc] 

Hhi 

[cm“^] 

IlHz 

[cm“^] 

0mol 

0CNM 

(Go = 1) 

<PCNM 
(Go = 5) 

140 

20+4 

17+5 

0.85 

2.7 

0.58 

80 

11+3 

49+5 

4.4 

1.5 

0.32 

40 

5+2 

101+8 

20.2 

0.7 

0.14 


Notes. Hi and H 2 particle density for different equivalent radii with 
corresponding Icr uncertainty, extracted using the information given in 
Fig.[T^and using equation[T^ The given uncertainties are the statistical 
uncertainties, but do not take into account the systematical uncertainties 
due to the saturation of the optical depth spectra. The model parameters 
and 0CNM are calculated using equations \n\ and|15[ respectively. 
(/irNM is calculated for solar values (Go = l,Z = 1) and for more realistic 
values of W43 (Go = 5, Z = 1.4). 


10^ 

I 

o 

a 



10 ^ 



biO 

O 


IQi 10^ 10^ 

^HI + ^H2 [MoPC^^] 


[Mac Low & Glover|2012[ ). 

To describe this transition, Krumholz et al. (2008 2009| l suggest 
an analytic model, which is based on the assumption of a spheri¬ 
cal cloud that is embedded in an isotropic external radiation field. 
They approximate that the transition between H i in the envelope 
and H 2 in the center occurs in an infinitely thin shell. An im¬ 
portant results is that an Hi column density of ~ IOMqPC^^ is 
necessary to shield the molecular hydrogen from the interstellar 
radiation field. If the cloud reaches this critical H i column den¬ 
sity the formation of molecular hydrogen is efficient enough that 
most of the hydrogen goes into molecular form and the column 
density of Hi remains relatively constant at this level. There¬ 
fore we should not expect to observe H i column densities larger 
than ~ lOMopc^^, contradicting our results presented in Fig. 
[T3| Furthermore, they show, that this H i column density thresh¬ 
old is independent of the external radiation field, but has a weak 
dependence on the metallicity of the gas. 

This model has three free parameters. First, the ratio of the mea¬ 
sured cold neutral medium (CNM) density hcnm to the minimal 
CNM density Wmin; 


«CNM - 0CNM ■ 


(15) 


The minimal CNM density is determined by the pressure balance 
with the warm neutral medium (WNM) and has a typical value 
of Umin ~ 7 - 8cm“^. As the range of pressure balance between 
the CNM and the WNM is limited, the maximum CNM density 
can be at most ~ 10 ■ Hence, 0 cnm can vary between 1 and 
10. |Kruniholz et aL] ( |2009| l assume 0 cnm ~ 3 for their fiducial 
value. The second free parameter is the ratio of the CNM density 
«CNM to the molecular density 


Umol — 0mol ' UCNM- 


(16) 


Krumholz et al. ( |2009) ) suggests that this ratio is (p„ioi ~ 10 


and should not vary significantly. The last free parameter is the 
metallicity Z. 

This model is supported by extragalactic observations 
( [Krumholz et al.||2009|l as well as observations in the Perseus 
Molecular Cloud by Lee et aL] ( |2012| ). They find a uniform 
Hi column density of Nhi ~ 6 - 8Mopc“^ for H 2 column 
densities up to Nh^ ~ 80Mopc“^. To fit their data to the model 
of jKrumholz et al.] ( |2009| l, they fixed i^nioi = 10 and Z - 1.0 Zq 
and fitted 0 cnm which reveals values of 0 cnm ~ 6 - 10. 

W43 is a more extreme test case for this theory as we have a 
much larger mass reservoir as well as an OB cluster that acts 
as a strong radiation source in the center. For this analysis, we 
extracted the column density of H i and H 2 for each pixel. As 


Fig. 17. The column density of the atomic hydrogen (visj = 60- 
120 km s“‘) as a function of the total hydrogen column density (Nm 
+ Nhj). The column density for the envelope (r > 140pc, Ah, = 
107 ± 20 Mo pc“^, Ah 2 = 116 ± 18 Mo pc“^) is subtracted. The solid and 
dashed lines describe the theoretical model of jKrumholz et alT]j2008| 
[2009) for solar values (Go = L Z = 1) and for more realistic values of 
W43 (Go = 5, Z = 1.4), respectively. The blue and black colors rep¬ 
resent the model parameters at different equivalent radii of r = 40 and 
140 pc, respectively. The black dotted line shows the fitted curve of the 
elliptical annuli analysis introduced in Sect.|5.3| 


a basis for the Hi column density we chose the optical depth 
corrected version with a velocity range of visr = 60- 120km s“'. 
For the H 2 column density, we again use the HiGal data (see 
Sect. |2.4| l. Even though the data base is the same as that 
used in Sect. HU we stress that the method is different. Here 
we conduct a pixel by pixel comparison, whereas in Sect. 
[5.5 1 we averaged the column density over elliptical annuli. 
In addition, we will focus merely on the inner part of W43 
with r < 140 pc and therefore, we have to subtract the col¬ 
umn density of the envelope. We assess the column density 
at r = 140 pc as the envelope column density and subtract 
this value. The column density for the Hi and H 2 envelopes 
are Ah, = 107 + 20Mopc“^ = 1.3 + 0.3 x lO^^cm^^ and 
Ah, = 116 + 18Mopc“^ = 7.3 + 1.1 x 10^^ cm“^, respectively. 
Figure shows the result of the pixel by pixel comparison of 
the HI and H 2 column density. Since we subtracted the column 
density of the envelope for this analysis, the values are smaller 
than in Fig. For better readability, we do not show each 
single pixel comparison, but we performed a pixel binning. We 
do not observe the predicted threshold for the H i column density 
of ~ 10 Mq pc“^, instead our data show high H i column density 
values, peaking between Ah, = 50 - 80Mopc“^. Nevertheless 
we will try to fit the analytic model of |Krumholz et al. ( 2009| l to 
our data. 

As described in Sec. HU we use elliptical annuli to estimate 
the particle density. We can use this information (Fig. 14 1 to 
determine the model parameters 0 cnm and ^moi using equations 
[^and[^ respectively. Given that the Hi and H 2 densities vary 
as a function of equivalent radius, the model parameters <^cnm 
and 0moi will also vary. We calculate the model parameters for 
three different distances r - 40, 80 and 140 pc that characterize 
the molecular dominated interior, the transition region and the 
well mixed outer area. To calculate <^cnm we have to know timin 
which is given by Krumholz et al. (20091 as: 


31 


1 -h3.1Z“« 


cm 


-3 


(17) 
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where Go and Z are the far-ultraviolet radiation intensity (in units 
of the Habing field) and the total metallicity, both normalized 
to their values in the solar neighborhood. For W43, we have a 
slightly higher metallicity of Z = 1.4Zo (|Motte et al.|2014 1 and 


a large radiation field which could be up to > 100 Go (Beuther 
et al. [20141) in the central region. We assume a more moderate 


radiation field for the outer regions of W43 with values around 
5 - IOGq. In Table|^ we present the Hi and H 2 density as well 
as the model parameters for different distances. The parameter 
0CNM is calculated for solar values (Go = 1, Z = 1) and for 
more realistic values of W43 (Go = 5, Z = 1.4) that results in 
Hmin = 7.5cm“^ and - 34.4cm“^, respectively. The solid 
blue and black lines in Fig.[^show the theoretical column den¬ 
sity for solar values and for the equivalent radius of r = 40 and 
140 pc, respectively, whereas the dashed lines represent the more 
realistic values of W43 (Go = 5, Z = 1.4). The dashed black 
line, which represents the outer area of W43 (r = 140 pc) with 
a moderate radiation field (Go = 5, Z - 1.4) might fit the data. 
However, the model parameters might have unrealistically low 
values of 0 cnm = 0.58 and 0moi = 0.85. Using a stronger radia¬ 
tion field will increase n^in and therefore decrease 0 cnm to even 
lower values. Hence, values with smaller distances and/or higher 
radiation fields predict column densities that are too high. 

We do not have a conclusive answer why the model suggested 
by lKrumholz et al!] ( |2008[[2009| l does not describe W43, but we 
suggest that the strong UV radiation produced by the central OB 
cluster ( |Blum et al.|1999| l and several further clusters in the en¬ 
vironment are responsible for the dissociation of the molecular 
hydrogen. Another explanation was given by Motte et al. (2014 1 , 
They performed a similar analysis and their conclusion was that 
the analytical model by |Krumholz et al.] ( |2008|[2009[ ) cannot de¬ 
scribe a complicated molecular cloud complex, as we see several 
transition layers between H i and H 2 along a single line of sight 
and the assumption of a simple spherical cloud, without internal 
radiation sources breaks down. 


6. Conclusions 

The Hi, OH, Recombination Line Survey of the Milky Way 
(THOR) is a Galactic plane survey covering a large portion of 
the Galactic disk (/ = 15 - 67°, \b\ < 1°). We use the VLA to ob¬ 
serve the 21cm Hi line, 4 OH lines, 19 Ha recombination lines 
and the continuum from 1-2 GHz. In this paper we present the H i 
data of the pilot field centered on the GMC associated with W43 
(/ = 29.2 - 31.5°, \b\ < 1°). The main results can be summarized 
as: 

1. We measured the average spin temperature of the neutral hy¬ 
drogen gas along the line of sight towards eight strong con¬ 
tinuum sources. Half of them are Galactic and half of them 
are extragalactic. We find a median value of Ts = 95.7 K, 
which is in agreement with other studies. 

2. We can estimate the optical depth for the Hi line toward 
strong continuum sources at various locations in W43. The 
measured optical depth saturates for the main velocity com¬ 
ponent of W43 (visr = 80 - lOOkms^*) with lower limits of 
T ~ 2.7 in the center. Hence, the derived Hi masses based 
on optically thin emission strongly underestimates the hy¬ 
drogen content. Employing further corrections for the weak 
and diffuse continuum emission, we obtain a lower limit for 
the HI mass of M - 6.6_i.8 x 10^ Mq for a velocity range of 
visr = 60 - 120km s“' and an area of / = 29.0 - 31.5°and 
\b\ < 1°. This is a factor of ~2.4 larger than the Hi mass 
estimates with the assumption of optically thin emission. 
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3. The measured Hi column density exceeds values of Ahi ~ 
150Mopc~^ over much of the inner region with r < 80pc. 
This is an order of magnitude larger than for low mass star 
forming regions such as Perseus. 

4. As we corrected the Hi column density for optical depth 
effects and the weak continuum emission, we are able to 
study the H i distribution spatially even in the innermost part 
of W43. We assumed an elliptical layered structure for the 
GMC associated with W43. This allows us to reconstruct the 
particle density of H i and we find a linearly decreasing H i 
density toward the center of W43 with values from hhi = 20 
to almost 0cm~^. Furthermore, we compared our results to 
the density of molecular hydrogen based on Herschel con¬ 
tinuum data. The density of molecular hydrogen shows an 
exponential increase toward the center of W43 with values 
rising from mhi = 15 to 200 cm~^. For smaller clumps, the 
density of H 2 can rise to even higher values. 

5. We compared our measurements to the analytic model sug¬ 
gested by |Krumholz et ak] P008| 2009| l. Our data does not 
show a sharp transition between H i and H 2 
find the predicted threshold for the H 


neither do we 
column density of 


~ 10 Mq pc“^. Based on these models, larger H i column den¬ 
sities should not exist, as molecular hydrogen should form 
for such high H i column densities. To fit the model, we have 
to assume low model parameters, which may indicate that 
the model is not applicable in a region with such a high radi¬ 
ation field. We suggest that the addition of an internal radia¬ 
tion field from a central cluster may be required to describe 
the observations. Thus, this work has shown that the H i con¬ 
tent of W43 and its relation to H 2 challenges current models 
of H 2 formation. 
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